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Abstract —An array low-density parity-check (LDPC) code is a 
quasi-cyclic LDPC code specified by two integers q and m, where 
q is an odd prime and m < q. The exact minimum distance, for 
small q and m, has been calculated, and tight upper bounds 
on it for m < 7 have been derived. In this work, we study the 
minimum distance of the spatially-coupled version of these codes. 
In particular, several tight upper bounds on the optimal minimum 
distance for coupling length at least two and m = 3,4, 5, that 
are Independent of q and that are valid for all values of g > go 
where go depends on m, are presented. Furthermore, we show by 
exhaustive search that by carefully selecting the edge spreading 
or unwrapping procedure, the minimum distance (when g is not 
very large) can be significantly Increased, especially for m = 5. 


we present the results of an exhaustive search over all cutting 
vectors, showing that a careful selection can increase the 
minimum distance significantly, especially for m = 5. 

We remark that the iterative decoding performance of these 
codes, for instance, on additive white Gaussian noise channels, 
is typically dominated by minimum absorbing sets and their 
multiplicities. However, improved decoding methods and ulti¬ 
mately maximum-likelihood decoding will not be trapped in 
absorbing sets and thus overcome the shortcomings of standard 
iterative decoding, in which case the minimum distance (and 
its multiplicity) will become an important performance metric. 


I. Introduction 

In this paper, we consider array-based spatially-coupled 
low-density parity-check (SC-LDPC) codes as introduced in 
[1] and subsequently studied in [2] from the perspective of 
absorbing sets. Array LDPC codes were originally proposed 
by Fan in [3] and are specified by two integers q and m, where 
q is an odd prime and m < q. 

Spatial coupling was first introduced in the coding theory 
literature by Felstrom and Zigangirov in [4], where they pro¬ 
posed convolutional LDPC or SC-LDPC codes. These codes 
have very good belief propagation (BP) decoding thresholds 
due to the phenomenon of threshold saturation. For instance, 
in [5], it was shown that for binary memoryless channels, the 
BP decoding threshold saturates to the maximum a posteriori 
decoding threshold of the underlying ensemble. 

Since the original work by Fan [3], several authors have 
considered the structural properties (including the minimum 
distance c?min) of array LDPC codes (see, e.g., [6-11]). For 
high rate and moderate length, these codes perform well under 
iterative decoding, and they are also well-suited for practical 
implementation due to their regular structure [12,13]. 

In this work, we present several tight upper bounds on the 
optimal minimum distance of array-based SC-LDPC codes 
for coupling length at least two and m = 3,4, 5, that are 
independent of q and that are valid for all values of g > go 
where go depends on m. Also, these bounds (and intermediate 
results in their proofs) can in some cases be used to quickly 
remove large parts of the search space when searching for 
optimal cutting vectors. For small values of g (m = 3,4, 5) 
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IF Array LDPC Codes 


The array LDPC code C{q,m), with parameters g and m, 
has length g^ and can be defined by the parity-check matrix 
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Since the number of ones in each row of the matrix in (1) 
is g and the number of ones in each column is m, the array 
LDPC codes are (m, g)-regular codes. Furthermore, its rank 
is qm — m + 1, from which it follows that the dimension of 
C(g, m) is g^ — qm + m — 1. 

In [7], a new representation for H{q,m) was introduced. 
In particular, since each column of the parity-check matrix 
H{q,m) has m blocks and each block is a permutation of 
(1, 0,0,..., 0,0)^, where (•)^ denotes the transpose of its 
argument, we can represent each column as a vector of integers 
between 0 and g — 1, where 







i.e., the 1-positions are associated with the integers modulo q. 
Furthermore, it follows from (1) and the integer representation 
in (2) that any column in an array LDPC code parity-check 
matrix is of the form 

{x,x + y,x + 2y,... ,x + {m — l)t/)^ mod q (3) 

where {x,y) € [g]^ and, for notational convenience, for any 
positive integer L, [L] = {0, 1,... ,L — 1}. Thus, a column 
can be specified by two integers x and y. Also, note that since 
there are q^ distinct columns in an array LDPC code parity- 
check matrix, any pair {x,y) G [g]^ specifies a valid column. 

In the following, it is also convenient to consider the parity- 
check matrix in (1) as an m x g array of g x g permutation 
matrices with row group indices i G [m] and column group 
indices j G [g], from which it follows that each column can 
be uniquely indexed by the pair (j, k), where j G [g] is the 
index of the column group and k G [g] is the index within a 
column group. Likewise, each row can be uniquely indexed 
by the pair {i, k), where i G \m\ is the index of the row group 
and k G [g] is the index within a row group. 

III. Array-Based SC-LDPC Codes 

Array-based SC-LDPC codes can be constructed from ar¬ 
ray LDPC codes by a special type of edge spreading or 
unwrapping procedure [4,14] specified by a “cutting” vector 

C = (Co, • ■ •, Cm-i), where 0 < Co < Ci < ''' < Cm-i < g- 

The purpose of the cutting vector is to generate two parity- 
check matrices Hq and Hi, each of size mg x g^ and initially 
filled with zeros, as follows: 

• For each Ci, i € [m], the g x g permutation matrices in 
row group i and column group j, j < Ci, of H (g, m) are 
copied into the corresponding positions in Hq. 

• Similarly, for each (i, i G [m], the q x q permutation 
matrices in row group i and column group j, Q<j< q, 
of H{q,m) are copied into the corresponding positions 
in Hi. 

Example 1: If g = 5, m = 3, and ^ = (1, 2,4), the matrices 
Hq and Hi become 

■j 0 0 0 o' 

Hq= I P 0 0 0 

J p2 p4 p 0 


As for the uncoupled case, it is convenient to consider the 
convolutional parity-check matrix in (4) as a (L -f l)m x Lq 
array of g x g permutation and all-zero matrices in which mg 
consecutive rows (resp. g^ consecutive columns) are referred 
to as a row (resp. column) section. Within each row (resp. 
column) section, the rows (resp. columns) correspond to a 
row (resp. column) group. In summary, each column can be 
indexed by a triple {l,j, k), where I G \L] is the column section 
index, j G [g] is the column group index within column section 
I, and k G [g] is the column index within column group j in 
column section 1. Similarly, each row can be indexed by a 
triple {l,i,k), where I G [L -f 1] is the row section index, 
i G [m] is the row group index within row section /, and 
k G [g] is the row index within row group i in row section 1. 

Note that as for the uncoupled matrix in (1), all columns in 
H{q,m,L,C) contain exactly m ones and each row in row 
sections 1, 2,..., L — 1 contains exactly g ones. However, the 
rows in row sections 0 and L can have fewer than g ones. 

In the following, we will denote the binary linear code 
defined by the coupled parity-check matrix H{q,m, L,<^) in 
(4) by C{q,m,L,C) and its corresponding minimum (resp. 
stopping) distance by d{q, m, L, C) (resp. h{q, m, L, C)). 

It was shown in [2, Theorem 5] that the minimum distance 
d{q, m, L, C) of C(g, m, L, C) is bounded below by the mini¬ 
mum distance d{q,m) of C{q,m) for all cutting vectors 

Now, let us define the optimal minimum distance of an 
array-based SC-LDPC code as the highest minimum distance 
over all possible cutting vectors as the coupling length L tends 
to infinity, i.e., 

rfopt(g,m)= lim dopt(g,TO,L) 

L—^oc 

where 

doptiq,m,L) = maxd{q,m, L,C). 

IV. Upper Bounds on dopt(g, m , L) 

In this section, we derive upper bounds on (iopt(g, w, L) 
for L > 2 and m = 3,4, 5, that are independent of g and that 
hold for all values of g > go where go depends on m. 

Now, define the range of a sorted (in nondecreasing order) 
sequence of column group indices ... ,g®Li) of 

length n as 
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Furthermore, define 
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For a given positive integer L (the coupling length) and a 
cutting vector an array-based SC-LDPC code is defined by 
the parity-check matrix 

Ho 

Hi Hq 

H{q,m,L,C)= Hi ■■. (4) 

Hq 

Hi_ 

of size {L + l)mq x Lq^. 


e{C) = max ((Co-ri) mod m - Q) mod g) 

jG[m-l] 

Jmax(C) = argmax ((Co-ri) mod m - Cj) mod g) 
je[m-l] 

for a given cutting vector 

Example 2: Let n = 4, g = 7, and jS*' = (0,1, 2, 5). Then, 
the range of jS’' is 

min (0 - (1 - 7 • 1), 1 - (2 - 7 • 1), 2 - (5 - 7 • 1), 

5 - (0 - 7 • 0)) -f 1 = min(6,6,4, 5) -f 1 = 5. 



A. The Case to = 3 


Lemma 1: The minimum distance dopt( 9 , 3,L), for L>2 
and q > 13, is upper-bounded by 6 . 

Proof: The proof is based on the template codeword from 
[7, Theorem 4]. For convenience of the reader, we restate the 
corresponding template support matrix below 


0 0 2i—2k 2i—2k —2i —2i 

0 —2i-\-k 0 —i —i —2i-\-k 

0 —4i-\-2k —2i-\-2k —42+2A; 0 —2i-\-‘2k 


(5) 


where q > 5, i G [g] \ {0}, and k G [g] with k f i, 2i, 
and where all operations are taken modulo q. The template 
support matrix is obtained by extracting the columns of 
the parity-check matrix corresponding to the support set of 
the underlying codeword. When constructing this matrix, the 
integer representation of the columns from (3) is used. 

Now, a column of the general form {x,x + y,x + 2y)’^ mod 
q, x,y G [g], in this matrix has column group index of y. 
Thus, the sequence of column group indices corresponding to 
the template matrix in (5) is 


(0, —2i -f k, —2i + 2k, —3i + 2k, i, k) mod q. 


For i = 1 and /c = 0, we get the following sorted sequence 


(-3,-2,-2, 0,0,1) mod g. ( 6 ) 


Since array LDPC codes are quasi-cyclic with period g, 
cyclically incrementing a sequence of column group indices 
jS'' = of length n of a codeword in an array 

LDPC code by an integer k > 1 results in the valid sequence 
Uo’^ + K, ... ,jn’Li + k) mod g of column group indices of 
a codeword that is obtained by cyclically shifting the given 
codeword by nq positions to the right. 

If the range of a sorted (in nondecreasing order) sequence 
of column group indices corresponding to a codeword in 
the uncoupled array LDPC code is at most e(C), it can 
always be cyclically incremented so that the corresponding 
codeword has coordinates with column group index j within 
[Cj„ax(C)>Cjm,,,(C)-ti) when < to - 1 = 2 , and 

within [Cm-i,g) U [0,Co) when jmax(C) = to - 1 = 2 . 
Consequently, there will exist a codeword in the coupled code 
(for all L > 2) with all coordinates within the same column 
section when jmax(C) < to— 1 = 2 , or within two consecutive 
column sections when jmax(C) = to — 1 = 2 . Thus, since 
min^e(C) = \^~\ = [§], there will exist a codeword in the 
SC-LDPC code of weight 6 for all cutting vectors as long 
as 

> range ((—3, —2, —2, 0,0,1) mod g). 

Since range((—3, —2, —2,0, 0,1) mod g) = 5 for g > 7 and 4 
for g = 5, the smallest g that satisfies this inequality is g = 13, 
and the result follows since the range of ( 6 ) for g > 13 is 
constant. ■ 

Example 3: Consider the case of g = 7, to = 3, and the 
template support matrix in (5) for i = 1 and k — 0. The 
sorted (in nondecreasing order) sequence of column group 
indices is (0,0,1,4, 5, 5) which has range 5. By cyclically 



incrementing this sequence by two, we get (0, 0, 2, 2, 3, 6 ). The 
corresponding support matrices (both of which correspond to 
codewords) are 

'0022551 , [ 002255 ' 

0 5 0 6 6 5 and 202110 
035305 J [402042 


respectively. For the cutting vector ^ = (4,5, 6 ), we have 
Jmax(C) = TO - 1 = 2, e(C) = 5, [Cm-I, g) u [0, Co) = [ 6 , 7) U 
[0,4) = {0,1, 2, 3, 6 }, and 
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From the highlighted blue columns in the matrix above we can 
identify a weight -6 codeword corresponding to the sequence 
(0,0, 2, 2,3, 6 ) of column group indices mentioned above. For 
other cutting vectors like C = (3,4, 6 ), we get e(C) = 4, which 
is less than the range, and there will be no codewords of the 
type in (5) with i = 1 and k = 0 (there are for i = 2 and 
fc = 0 ) in the coupled code. 

By combining Lemma 1 with [2, Corollary 6 ], we get the 
following theorem. 

Theorem 1: For L >2 and g > 13, the minimum distance 
t^opt(g 5 3, L) = 6 . 

B. The Case to = 4 

Lemma 2: The minimum distance (iopt(g, 4, L), for L >2 
and g > 41, is upper-bounded by 10. 

Proof: The proof follows the same main idea of the proof 
of Lemma 1 above using the template support matrix (g > 11) 
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which we have found using the algorithm from [10]. The 
remaining technical details are omitted for brevity. ■ 

Note that the template support matrix for to = 4 given in 
Fig. 3 in [ 8 ] will fail to prove the above result, since the sorted 
(in nondecreasing order) sequence of column group indices has 
range 19 for g > 23. As a consequence, the template matrix in 
Fig. 3 in [ 8 ] can only prove an upper bound of 10 for g > 73. 

Proposition 1: There is a codeword of weight 10 in 
C(g,4,L,C), g 13, for all positive integers L > 2 and 
cutting vectors ^ if one of the following conditions are true. 

1) Cl ~ Co > 2 and C 2 — Ci ^ 9- 

2) Cl - Co > 9 and C 2 - Ci > 2. 

3) C 3 - C 2 > 2 and g + Co - Ca > 9. 

4) Ca - C 2 > 9 and g + Co - Ca > 2. 

5) max(Ci — Co,C 2 — Ci,Ca — C 2 ,g + Co — Ca) > H- 

Proof: The proof is omitted due to lack of space. ■ 

Lemma 3: The minimum distances (iopt(31,4, L) and 
(iopt(37,4, L), for L > 2, are both upper-bounded by 10. 

Proof: For the case g = 31, the number of cutting vectors 
that satisfy none of the five conditions in Proposition 1 is only 




















35, while the total number of possible cutting vectors is 35 960. 
Thus, by running the algorithm from [15], adapted to the case 
of SC codes, on these 35 cases we have verified that there are 
indeed codewords of weight 10 for all cutting vectors, and the 
result follows. 

For q = 37, there are no cutting vectors that satisfy none 
of the five conditions in Proposition 1, from which it follows 
that fiopt(37,4, L) < 10 for all positive integers L>2. ■ 

By combining Lemmas 2 and 3, [2, Theorem 5], and [8, 
Corollary 4.2], we get the following theorem. 

Theorem 2: For L > 2 and q > 31, the minimum distance 
t^opt 4, Z/) = 10. 

C. The Case m = 5 

Theorem 3: The minimum distance (iopt (95 5, L), for L >2 
and q > 59, is upper-bounded by 12. 

Proof: The proof follows the same main idea of the proofs 
of Lemmas 1 and 2 above using the template support matrix 
from Fig. 4 in [8], and is omitted for brevity. ■ 

Note that there appears to be no equivalent to Proposition 1 
(except for the last condition) for m = 5, since the template 
support matrix from Fig. 4 in [8] (which is used for the proof 
of Theorem 3) does not have the required structure. Also, the 
algorithm from [10] was not able to identify other nonequiv¬ 
alent template support matrices with a suitable structure. This 
has also been “confirmed” by the fact that for g = 53 we have 
been able to identify a cutting vector that gives a minimum 
distance of 14 (see Table III). 

Theorem 4: The minimum distance dopt{q, 5, L), for L >2 
and q > 29, is upper-bounded by 16. 

Proof: The proof follows the same main idea of the proofs 
of Lemmas 1 and 2 and Theorem 3 above using a template 
support matrix which was found using the algorithm from [10]. 
Due to lack of space, the actual template support matrix and 
the technical details of the proof are omitted. ■ 

V. Numerical Results 

In this section, we present some numerical results for the 
cases TO = 3, 4, and 5. 

A. The Case to = 3 

By performing an exhaustive search over all cutting vectors 

the optimal minimum distance (iopt(9j3) was determined 
for 5 < g < 11 (see Table III). For each of these values of 
q, the dniin-optimal cutting vector is not unique, i.e., there are 
several cutting vectors that give the best minimum distance. 
For instance, for g = 5, we found (for L = 10) the hmin- 
optimal (where hmin denotes stopping distance) cutting vectors 
displayed in Table 1. For each cutting vector, in the second 
and fourth columns we have tabulated the corresponding min¬ 
imum and stopping distance, respectively. The corresponding 
multiplicities are tabulated in the third (minimum distance) 
and fifth (stopping distance) columns. Note that the optimal 
cutting vector ^ = (2,3, 5) for (3, 3) and (4, 2) absorbing sets 
from [2] is also optimal when it comes to minimum/stopping 
distance, since it gives the optimal minimum/stopping distance 


TABLE I 


Minimum/Stopping Distance Results for Array-Based SC-LDPC 

Codes for q = 

5, m = 3, L = 

10, AND Different hmin-OPTiMAL 


Cutting Vectors 


C 

d(5.3. 10, C) 

Mult. 

?i(5.3.10,C) 

Mult. 

(0.1.3) 

10 

20 

10 

65 

(0.2.3) 

10 

20 

10 

65 

(1.2.4) 

10 

19 

10 

59 

(1.3.4) 

10 

19 

10 

59 

(2.3.5) 

10 

20 

10 

65 

(2.4.5) 

10 

20 

10 

65 


TABLE II 



Minimum/Stopping Distance Results for Array-Based SC-LDPC 

Codes for q = 

7, m = 4, L = 

10, AND Different /imin-OPTiMAL 


Cutting Vectors 


C 

<i(7,4,10,C) 

Mult. 

d(7,4,10,c) 

Mult. 

(0. 2.3.5) 

14 

30 

14 

401 

(0.2.4. 6) 

14 

30 

14 

695 

(1.3.4.6) 

14 

29 

14 

393 

(1.3. 5. 7) 

14 

30 

14 

695 

(2.4. 5.7) 

14 

30 

14 

401 


of 10. However, it gives slightly more stopping sets of size 
10 (see Table I). On the other hand, for q = 7, the optimal 
cutting vector C = (2,4,6) for (3,3) and (4, 2) absorbing sets 
from [2] is not optimal when it comes to minimum distance, 
since it gives a minimum distance of only 6 (the optimal 
minimum distance is 8). For q = 5,7, dmin-optimal cutting 
vectors that also give the lowest possible minimum distance 
multiplicity, denoted as ((inii„, -optimal cutting vectors, 

are displayed within the parentheses in the sixth column of 
Table III (first and second row, respectively). 

For <7 > 13, it follows from Theorem 1 that the optimal min¬ 
imum distance is 6. Also, all cutting vectors give a minimum 
distance of 6, but possibly with different multiplicities. 

B. The Case to = 4 

By performing an exhaustive search over all cutting vectors 
^ (with the help of Proposition 1 to reduce the search space 
when q > 11) the optimal minimum distance dopt(( 7 ,4) was 
determined for 5 < q < 29. For each of these values of q, 
the dmin-optimal cutting vector is not unique, i.e., there are 
several cutting vectors that give the best minimum distance. 
For instance, for q = 7, we found (for L = 10) the hmin- 
optimal cutting vectors displayed in Table 11. For each cutting 
vector, in the second and fourth columns we have tabulated the 
corresponding minimum and stopping distance, respectively. 
The corresponding multiplicities are tabulated in the third 
(minimum distance) and fifth (stopping distance) columns. For 
5 < <7 < 29, the optimal minimum distance and corresponding 
(rfmiii) -optimal cutting vectors (displayed within the 

parentheses) are given in the fourth column of Table III (the 
first eight rows, respectively). For q > 31, it follows from 
Theorem 2 that doptiqA) = 10. 

C. The Case to = 5 

For <7 = 7, we have performed an exhaustive search over 
all possible cutting vectors. The optimal minimum distance 
rfopt(7, 5) is 14 and a (dmin, -optimal cutting vector 

is displayed within the parentheses in the second column of 



TABLE III 

Optimal Minimum Distance Results eor Array-Based SC-LDPC Codes eor Dieeerent Values oe q and m. Values in Bold Are New 
Results, While Non-Bold Values Are Taken From the Literature. The Results Are Compared With Those oe Array-Based 
Uncoupled LDPC Codes Taken From [7, 8,10], Within the Parentheses the Corresponding Cutting Vectors Are Displayed. Cutting 
Vectors in Italics Are dmin-OPTiMAL, While Cutting Vectors in Bold Are (dmin, Afd^.^)-OPTiMAL 


Q 

dopt(9) 6) 

d(q, 5) 

dopt (<?? 4) 

d(q,A) 

dopt3) 

d(q, 3) 

5 



10 (1, 2, 3, 4) 

8 

10 (1, 2, 4) 

6 

7 

14 (1,2, 3,4. 6) 

12 

14 (1,3, 4, 6) 

8 

8 (1, 2, 5) 

6 

11 

22 (0,1,2, 5, 8) 

10 

14 (1, 3, 6, 8) 

10 

6 

6 

13 

24 (0, 2, 5, 7, 10) 

12 

12 (0, 3, 7, 9) 

10 

6 

6 

17 

20, 22, or 24 (0,4,7, 11, 14) 

12 

12 (2, 5, 10, 13) 

10 

6 

6 

19 

20 or 22 (0,4, 8, 12, 16) 

12 

12 (3, 7, 12, 16) 

10 

6 

6 

23 

20 (5, 10, 15, 19, 23) 

12 

12 (4, 8, 15. 19) 

10 

6 

6 

29 

16 (8, 11, 18, 26, 29) 

12 

12 (4, 11, 18, 25) 

10 

6 

6 

31 

16 (9, 15, 21,24,31) 

12 

10 

10 

6 

6 

37 

16 (0, 10, 19, 28, 30) 

12 

10 

10 

6 

6 

41 

16 (0, 9, 18, 25, 32) 

12 

10 

10 

6 

6 

43 

16 (0, 8, 19, 24,33) 

12 

10 

10 

6 

6 

47 

16 (7, 16, 26, 36,47) 

12 

10 

10 

6 

6 

53 

14 (0, 11,22, 33, 42) 

12 

10 

10 

6 

6 

59 - 79 

12 

12 

10 

10 

6 

6 

> 79 

10 or 12 

10 or 12 

10 

10 

6 

6 


Table III (second row). For q = 11 and 13, the optimal 
minimum distance is as high as 22 and 24, respectively, and 
rfmin-optimal cutting vectors are shown within the parentheses 
in the second column of Table III (third and fourth row, 
respectively). For g = 17 and 19, upper and lower bounds on 
the optimal minimum distance, tabulated in the second column 
of Table III, have been established from an exhaustive search 
over all cutting vectors. The displayed values are bounds since 
we were not able to exhaustively enumerate all codewords of 
weight at most r, when t > 20, for a given cutting vector, for 
these two values of q. Within the parentheses a corresponding 
cutting vector (which establishes the upper bound) is also 
tabulated. The lower bounds were determined by running the 
algorithm from [15], adapted to the case of SC codes, for the 
specific cutting vectors displayed within the parentheses. For 
23 < g < 53, the tabulated values (in the second column of 
Table III) are again the exact values of (iopt(ii, 5) and also 
(within the parentheses) dmin-optimal ((dmin, 7Vd,„j„)-optimal 
for q = 53) cutting vectors are displayed. For 59 < <7 < 79, 
do-ptiq, 5) = 12, which follows from Theorem 3, [10, Table I], 
and Theorem 5 in [2]. For g > 79, it follows from Theorem 3, 
Corollary 4.4 in [8], and Theorem 5 in [2] that (iopt(9, 5) is 
either 10 or 12, although we conjecture it to be 12. 

VI. Conclusion and Future Work 

In this work, we have studied in detail the minimum distance 
of array-based SC-LDPC codes. Several tight upper bounds 
on the optimal minimum distance for coupling length L > 
2 and m = 3,4, 5, that are independent of q and that are 
valid for all values of <7 > <70 where go depends on m, have 
been presented. Furthermore, we have conducted an exhaustive 
search over all cutting vectors for small values of q {m = 
3,4, 5) which shows that by carefully selecting the cutting 
vector, the minimum distance (when q is not very large) can 
be significantly increased, especially for m = 5. 

An interesting topic for future work is to consider the corre¬ 
lation with absorbing sets. In particular, will a (dmin, -^dmin)" 
optimal cutting vector also be close-to-optimal when it comes 
to problematic absorbing sets, and/or vice-versa? 
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